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ABSTRACT 


The  calculations  presented  herein  are  based  on  a  first -collision 
model  which  allows  collisions  between  the  free-stream  molecules  and 
the  molecules  which  are  re -emitted  from  the  surface  of  the  body.  The 
net  effect  of  these  collisions  is  to  partially  shield  the  body  from  the  free 
stream,  reducing  both  the  drag  and  heat-transfer  coefficients  from  the 
corresponding  values  experienced  in  free-molecule  flow.  Free-stream 
Mach  number  is  taken  to  be  essentially  infinite  and  the  molecules  are 
assumed  to  be  re-emitted  from  the  body  surface  at  the  most  probable 
velocity,  instead  of  possessing  a  velocity  distribution.  These  assump¬ 
tions  enable  the  distribution  of  flux  to  a  given  body  to  be  expressed  as  a 
function  of  the  degree  of  rarefaction,  as  represented  by  the  appropriate 
Knudsen  number.  This  Knudsen  number  is  composed  of  a  characteristic 
body  dimension  and  the  mean  free  path  of  the  re-emitted  molecules  rela¬ 
tive  to  free-stream  molecules.  The  integral  equation  expressing  the 
incident  flux  distribution  on  a  general  body  is  developed,  and  solutions 
are  presented  for  the  disk  normal  to  the  free  stream  and  for  sharp  cones 
of  various  apex  angles,  at  zero  angle  of  attack.  These  flux  distributions 
are  then  integrated  to  give  drag  coefficients  ratioed  to  the  corresponding 
free-molecule  values. 
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NOMENCLATURE 


A  Area 

man  Total  shielding  of  point  located  at  m  resulting  from  a  unit 

flux  density  re-emitted  from  a  concentric  ring  of  radius  r 

b  Nondimensional  distance  from  point  Si  to  point  S  in 

plane  of  disk 

C  Constant  in  Eq.  (1) 
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Denotes  particular  point  for  which  the  flux  density  is 
desired 
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re -emitted 

Point  defined  in  Section  2.  1 
Pressure 
Total  heating  rate 
Heat -transfer  rate  per  unit  area 
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Radial  location  of  point  S 
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SECTION  l 
INTRODUCTION 


This  paper  is  addressed  to  the  calculation  of  the  interaction  of  a 
gas  flow  with  a  body  in  the  transition  flow  regime  between  continuum 
and  free-molecule  flow.  Several  attempts  have  been  made  to  extend 
continuum-flow  calculations  into  the  transition  regime,  the  best  known 
of  which  are  perhaps  by  Van  Dyke  (Ref.  1)  and  Cheng  (Ref.  2).  Cheng 
has  solved  exactly  a  set  of  simplified  equations  which  approximately 
describe  the  flow  field.  The  model  used  becomes  less  appropriate  as 
higher  Knudsen  numbers  are  approached,  but  the  results  obtained 
approach  the  correct  limit  for  free-molecular  flow.  Van  Dyke's  second- 
order  theory  is  valid  for  some  portion  of  the  transitional  regime,  close 
to  continuum  flow,  but  becomes  progressively  more  inaccurate  as  the 
flow  becomes  more  rarefied.  In  principle  a  "third'1  or  higher  order 
theory  could  be  proposed  to  extend  the  region  of  validity  to  higher  Knudsen 
numbers,  but  the  formulation  of  these  higher  order  terms  becomes  very 
unwieldly  and  then  there  is  no  assurance  that  they  would  tend  to  converge. 

All  of  the  aerodynamic  forces  acting  on  a  body,  as  well  as  the  gross 
heating  effects,  are  given  by  the  summation  of  the  local  interactions  of  the 
flow  with  the  body  surface.  .  These  local  interactions  can  be  examined  from 
the  microscopic  point  of  view  as  local  collision  processes  between  the  wall 
and  the  molecules  of  gas  in  the  vicinity  of  the  wall.  If  the  velocity  distribu¬ 
tion  and  the  interaction  potential  between  the  molecules  are  known,  the 
local  flux  of  momentum  and  energy  to  the  wall  can  be  determined. 

The  problem  with  this  approach  is  that  the  state  of  the  gas  molecules 
in  the  vicinity  of  the  wall  is  not  known  and  cannot  easily  be  determined 
except  in  the  special  case  of  free-molecule  flow.  In  this  report  the  solu¬ 
tion  to  the  problem  is  formulated  in  such  a  manner  that  the  transition 
regime  is  entered  from  the  free-molecule  flow  regime. 

The  solution  presented  herein  is  based  on  a  first -collision  model, 
developed  in  Ref.  3,  which  allows  collisions  between  the  free-stream 
molecules  and  the  molecules  which  are  re-emitted  from  the  surface  of 
the  body.  After  a  collision,  the  participating  molecules  are  considered 
lost,  in  that  subsequent  events  are  disregarded,  so  the  net  effect  is  a  partial 
shielding  of  the  body  from  the  free-stream  molecules  by  the  molecules 
re -emitted  from  the  surface. 

The  effectiveness  of  this  shielding  is  a  function  of  the  degree  of  rare¬ 
faction  of  the  flow  which  is  expressed  in  terms  of  an  appropriate  Knudsen 
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number.  This  Knudsen  number  is  based  on  the  mean  free  path  of  the 
re-emitted  molecules  relative  to  collisions  with  the  free-stream  mole¬ 
cules.  This  single  mean  free  path  value  is  determined  by  assuming 
that  all  incident  molecules  are  re -emitted  from  the  surface  at  the  most 
probable  velocity  instead  of  possessing  a  velocity  distribution. 

Once  the  distribution  of  incident  free-stream  molecules  on  the  sur¬ 
face  of  a  body  is  determined,  the  pressure  distribution  and  heat -transfer 
distribution  can  be  found.  These  distributions  can  then  be  integrated  to 
yield  drag  coefficient  and  total  heat  flux  to  the  body  in  near-free-molecule 
flow. 

For  any  experimental  data  a  problem  arises  concerning  the  deter¬ 
mination  of  the  effective  Knudsen  number,  Kn,  which  evolves  naturally 
in  the  theory.  Reference  3  gives  the  result 

Kn  =  Cy'/2  M„/Re..  (1) 

where  Kn  is  the  Knudsen  number,  is  the  free-stream  Mach  number, 
Rem  is  the  free-stream  Reynolds  number,  and  C  is  a  irfree  constant" 
which  is  chosen  to  correlate  the  experimental  data  with  theory.  This 
approach  yields  a  curve  which  passes  through  the  data,  whose  shape  has 
theoretical  justification,  and  whose  value  has  an  experimental  basis. 

While  this  approach  is  not  completely  satisfactory  from  a  theoretical 
viewpoint,  its  justification  can  be  argued  on  the  basis  that  this  constant 
is  necessary  to  account  for  the  uncertain  force  field  which  exists  between 
the  actual  molecules.  For  example,  for  billiard  ball  molecules  this  con¬ 
stant  has  a  theoretical  value  of  C  =  1.  26  when  calculating  free-stream 
Knudsen  number,  Kn,,.  In  Ref.  3,  which  deals  with  hypersonic  spheres 
in  near-free -molecular  flow,  it  was  found  that  C  -  2.  446. 

The  results  presented  herein  are  given  in  terms  of  the  rarefaction 
parameter,  R/X,  where  R  is  a  characteristic  body  dimension  and  X  is  the 
mean  free  path  of  the  re -emitted  molecules  in  the  uniform  free  stream. 
These  results  are  strictly  applicable  only  for  small  values  of  R/X.  How¬ 
ever,  the  theory  gives  the  proper  limit  for  flow  with  large  R/X  in  the  case 
of  a  blunt  body. 


SECTION  II 
ANALYSIS 


2.1  GENERAL  AXISYMMETRIC  BODIES 

The  following  development  assumes  that  the  free  stream  is  composed 
of  molecules  having  a  common  direction  and  velocity,  which  is  equivalent 
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to  assuming  an  infinite  Mach  number.  Consider  the  general  body  immersed 
in  such  a  stream  as  sketched  below. 


For  a  streamwise  tube  of  unit  cross-sectional  area  terminating  at  the  body 
surface  centered  on  point  S]_,  the  difference  in  number  density  of  free- 
stream  molecules  entering  the  tube  and  the  number  density  of  molecules 
actually  surviving  to  strike  the  surface  of  the  model  is  given  by 

-  MS.)]  =  J%dx  (2) 

O 

where  and  p*  (Sj)  are  the  number  density  of  molecules  at  infinity  and 
the  number  density  striking  the  surface  of  the  model,  respectively,  and 
ip  is  the  collision  number  density  per  unit  time  at  P.  To  find  p,  consider 
the  contribution  of  the  events  at  dA  on  the  surface  to  the  events  at  the  point 
P.  Then  one  sees  that 


dO  =  u  ,  (S)  •  COS  (J  dA  • 

OC  •  1 


cos  (f> 
TT 


7 


The  significance  of  the  various  quantities  is  as  follows: 


(3) 


=  number  flux  impinging  on  dA  (Pi(S)  is  the  number 
density  striking  the  surface  at  location  S. ) 

cos  cr  dA  =  unit  frontal  area 

=  probability  that  a  molecule  is  re-emitted 
along  H  (diffuse  reflection) 

~~  =  decrease  in  density  at  P  caused  by  spherical 
spreading  of  re-emitted  molecules 

_  i/x 

e  =  probability  of  surviving  passage  to  P  without 
suffering  a  collision 

=  probability  of  collision  per  unit  length  at  P 
dp  —  collisions  per  unit  volume  per  unit  time 
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0  is  then  given  by  integration  over  the  portion  of  the  body  surface  which 
can  be  "seen'1  from  P.  Thus, 


area  of  body 


A-/  / 


COS  <7 
^  =  0 


-  i/A. 
e 


dx  dA 


Defining  y  =  then 


(4) 


y(Si)  = 


(5) 


where  S  denotes  the  location  on  the  body  surface  from  whence  the  mole¬ 
cules  are  emitted  and  Si  denotes  the  location  on  the  body  from  which  the 
x  vector  is  drawn,  that  is,  the  location  on  the  body  surface  which  is  partly 
shielded  from  the  free  stream  by  collisions  between  molecules  from  S  and 
the  free  stream. 


Equation  (5)  is  valid  for  any  body  configuration,  and  the  solution  gives 
the  flux  density  which  impinges  on  the  point  Si  on  the  body  surface.  The 
inner  integral  is  a  function  of  particular  body  geometry  and  the  location 
of  Si  and  S.  This  equation  is  a  nonhomogenous  Fredholm  linear  integral 
equation  of  the  second  kind  (Ref.  4).  The  solution  was  obtained  by  dividing 
the  forward-facing  surface  of  the  body  into  small  increments  of  area  such 
that  y(S)  over  each  small  area  could  be  taken  to  be  constant.  This  allows 
y(S)  to  be  moved  outside  the  integral  and  Eq.  (5)  can  be  rewritten  as 


y (m/M )  =  1  - 


A  A(n)  x 

— —  y(n/M) ^  cos  (7  ^ 

n  =  0 


cos  d>  —  l/\  .  . . 

-  e  dx  dA 


(6) 


cos  ^  =  0 

where  M  is  the  total  number  of  area  increments  used,  and  n  and  m  denote 
the  particular  area  increments.  If  man  is  defined  by 

A  Ain)  x  =  oo 

cos  (f>  —  l/k 


1  C  C  cos 

-  I  cos  cr  I  - 

AJ  J 


dx  dA 


cos  =  0 


then  the  equation  can  be  expressed  as  the  following  array: 
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(0/M)  =  1  - 

oao  >'(0/M)  - 

Qa  j  y  (1/M)  - 

y (2/M)  - 

0a3  ■ 

■  0aM  y(1) 

(1/M)  -  1  - 

^  y  (O  ’M'I  - 

1a)  y(l  M)  - 

La2  v (2/M)  - 

I—1 

•  .  a  y  (1) 
1  M  ' 

,(1)  .  1  -  M.0  ylO.'M)  -  ,.1/M)  -  M.2  v(2,M)  -  yUUI).  .  .  ^  yll) 

► 

(l  -  ,(«•»!)  +  y(l  M)  +  0«2  v(2,M)  -  (S  y (3 'Ml  ♦  ...  y(lt  .  1 

taQ  y(0/M)  +  (!  +  xaj)  y(l/M)  +  yi(2/M)  +  . ^  y{l)  =  1 


MaO  >(0/!V,)  +  Ma)  >(1'M)  . . O  +  m“m)  y(1)  “  1 

Thus  the  solution  is  given  by 


In  principle  the  solution  could  be  obtained  in  this  manner  for  any  body 
geometry,  but  practically  the  number  of  points  is  limited  by  the  size  of 
the  M  by  M  man  matrix  to  be  inverted. 

The  a1  s  actually  represent  the  influence  of  the  n  location  on  the  m 
location  on  the  body.  For  both  free -molecule  and  continuum  flow,  the  off- 
diagonal  elements  are  zero.  This  is  because  in  free -molecular  flow  only 
an  insignificant  number  of  re-emitted  molecules  suffer  a  collision  in  a 
finite  distance  from  the  point,  while  at  the  continuum  limit,  R /X— all 
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the  significant  interactions  occur  in  the  immediate  vicinity  of  the  point. 
Thus,  nan  goes  uniformly  from  zero  at  R/A — *-0  to  unity  at  R/A  — giving, 
for  any  body  geometry,  y— - 1.  0  for  free-molecule  flow  and  y — -0.  5  for 

R/A—®. 


It  is  obvious  that  the  largest  coefficient  will  come  from  the  considera¬ 
tion  of  the  effect  of  an  element  upon  itself.  That  is,  when  m  and  n  coincide. 
In  this  situation,  H — ►x  and  cos  $ — ►cos  a,  and  taking  cos  a  dA  as  the  incre¬ 
mental  normal  area,  dAn,  the  contribution  to  nam  becomes 

A  An  x -*  oo 


0 


in  which  the  integral  is  unbounded  for  x — *-0.  The  remedy  is  to  choose  a 
coordinate  system  for  dAn  such  that  dAjf*-0  at  the  point  in  question.  A 
cylindrical  coordinate  system  with  the  origin  at  the  point  in  question  satis¬ 
fies  this  requirement  and  gives  a  finite  result  for  nam.  For  a  more  com¬ 
plete  discussion,  the  problem  must  be  expressed  in  a  particular  coordinate 
system. 


2.2  FLAT  DISK  NORMAL  TO  FLOW 


The  coordinate  system  used  for  this  case  is  cylindrical,  with  the  axis 
shifted  to  the  point  Si  at  which  the  flux  density  is  desired  so  that  the  incre¬ 
mental  normal  area,  dAn,  approaches  zero  as  S  and  Si  coalesce.  Using 
the  local  nondimensional  radius,  ri/R  for  Si  and  r/R  for  S,  the  equation 
to  be  evaluated  is 


y(ri/R) 


1 


1 


ormal  Area 
y  (r/R) 


-1/ A 

e 


dx  dAn 


(8) 


with  the  symbols  defined  by  the  sketch  below. 


6 


AE  DC-T  ^  -67-26 


For  this  body  shape  the  following  relations  may  be  found; 

cos  dx  =  Ak 

d.\„  =  R2  b  db  d/3 

and  as  cos  # — >-0,  x — -0  and  k  =  Rb.  Using  these  relations  and  assigning 
the  appropriate  limits  of  integration,  the  problem  becomes  the  solution  of 


where 


(r/'R) 


+ 


b 


2  (r,  /  R)  b  cos  /9 


Now  considering  the  innermost  integral,  one  finds 


E, 


where  En  denotes  the  exponential  integral  of  order  n.  Using  this  result, 
Eq.  (9)  becomes 


y(r,/R) 


If  y(r/R)  is  approximately  constant  for  a  small  interval  in  b. 


then 


b-A!><b<b  +  ^ 


2 


'  y(r/iH|b  -f  jE.  [-n.(b  -  ^)] 


(10) 
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The  contribution  of  area  dAn  =  b  db  d/3  located  at  S  to  the  shielding 
of  point  Si,  is  given  by 


y(r/H)s  DA 


—  A  /3  y(r/R)s  ]E, 


which  defines  DA. 


Let  rj/R  =  Sj  =  0.  lm  and  r/R  =  S  =  0.  In,  where  n  and  m  are  integers 
which  vary  between  0  and  10.  (This  assumes  that  the  disk  is  divided  into 
ten  concentric  rings  on  each  of  which  y(s)  =  constant. )  The  coefficient 
man  can  then  be  defined  as 

a  =  I  DA  (12) 

m  n  '  ' 

with 


r  =  constant  =  O.ln 
Tj  =  constant  =  0.1m 

This  nomenclature  is  clarified  by  the  following  sketch: 


The  frontal  area  of  the  disk  is  covered  by  small  units  of  area  as 
shown,  which  are  defined  by  a  cylindrical  coordinate  system  with  the 
origin  located  at  Si.  The  corresponding  DA  is  calculated  and  its  con¬ 
tribution  is  credited  to  the  proper  man.  The  location  of  Si  is  then 
changed,  so  that  a  complete  set  of  man  coefficients  is  obtained,  with 
m  and  n  taking  values  from  zero  to  ten.  This  gives  an  11-  by  -11  matrix, 
the  inversion  of  which  yields  the  flux  density  impinging  on  the  body  at 
eleven  evenly  spaced  radii. 
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A  simple  extension  of  the  calculation  has  been  included  to  give  the 
influence  of  the  presence  of  the  disk  on  the  fluid  flow  in  its  plane  beyond 
the  edge  of  the  disk. 


2.3  SHARP-NOSED  CONES 

For  this  configuration,  a  cylindrical  coordinate  system  is  used,  with 
the  coordinate  system  centered  on  the  cone  axis.  This  arrangement  is 
suitable  for  all  cases  in  which  the  points  S  and  Si  do  not  coalesce.  When 
the  points  S  and  Si  do  coalesce,  the  origin  of  the  coordinate  system  is 
shifted  to  the  point  itself,  but  only  for  the  calculation  of  the  influence  of 
flow  interactions  in  the  immediate  vicinity  of  this  point.  That  is,  after 
this  particular  DA  is  calculated,  the  coordinate  system  origin  is  immedi¬ 
ately  shifted  back  to  the  cone  axis,  and  the  calculation  proceeds.  The 
pertinent  geometry  is  shown  below. 


Here,  Si  denotes  the  location  of  the  point  at  which  the  shielding  is 
desired  because  of  the  influence  of  the  area  dA  located  at  S.  The  angle 
< i>  is  the  angle  between  the  normal  to  the  surface  (at  S)  and  £ ,  the  path 
of  particles  emitted  from  dA  which  may  collide  at  P  with  free-stream 
particles.  The  quantity  XL  is  the  lower  limit  on  the  x-integration  (for 
x — ►XL,  Si  can  no  longer  Tisee"  S).  When  cos  0  =0,  x  =  XL. 
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Taking  the  cone  half-angle  as  6C,  the  necessary  quantities  are  as 
follows : 

=  [x  -  cos  dc  (d,  -  d)  1  +  sin3  dc  (d,3  +  d2  -  2dd,cos  <y ) 

XL  =  d,  cos  0C  {1  -  cos  Ci>  ) 

COS  </>  =  -Sln  [x  -  di  cos  8C  (1  -  COSO))]  ^ 

z 

where  u  is  the  angle  between  the  plane  through  the  cone  axis  containing 
Si  and  a  plane  containing  both  S  and  the  axis.  These  relations  in  Eq.  (13) 
are  used  to  yield 


Normal  Area 


-i/\ 

cos<£  — -  dx  dAn 

t 2 


(14) 


where 


dAn  =  d  sin  dc  d  0)  d(d  sin  6C) 


It  is  convenient  to  nondimensionalize  Eq.  (14)  by  the  cone  base  radius, 
R,  which  is  given  by  R  =  Smax  sin  0C.  Using  this  result  and  retaining  the 
previous  symbols  for  the  new  nondimens ional  quantities,  the  equation  for 
DA  becomes  (when  S  and  Sj  do  not  coalesce) 


DA(S  4  S,) 


_2 

n 


4  o  -  (d)  (Ad) 
X 


COS  0 


-(RAH 


-e 


dx 


(15) 


This  integral  can  be  evaluated  numerically  with  reasonable  accuracy  using 
a  computer,  provided  the  integration  is  performed  over  an  appropriate 
interval.  The  procedure  used  was  to  calculate  approximately  the  appropri¬ 
ate  interval,  integrate  the  function  through  this  interval  and  then  check  to 
see  if  further  integration  was  required.  Finally,  the  value  of  the  integral 
from  this  upper  limit  to  infinity  was  approximated,  using 

cos  (p  sin  6C  ,  dx  =  di 

and 

oo 

cos  - -  dx 

tx 

X 

upper 


CO 

/ 


■  (R/X)/ 


di  = 


oc 

sin  C 

upper  J 


-(R/A)z  fur 


dz 


upper 


sin  0 c 

^upper 


(16) 
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Thus, 


DAis  *  v  - 


upper 


(17) 


When  S  and  coalesce,  the  shielding  of  the  point  by  its  own  re- 
emitted  molecules  is  computed  by  assuming  the  element  of  normal  area, 
dAn,  is  actually  oriented  normal  to  the  free  stream.  This  unit  of  area 
is  subdivided  into  smaller  units  as  shown  in  the  sketch,  which  shows  an 
axial  view  of  the  cone.  A  cylindrical  coordinate  system  is  used  with  the 


origin  at  the  points  S  and  S]..  The  angle  £  is  varied  between  0  and  180  deg 
by  equal  increments  of  A/3,  the  appropriate  z  is  calculated,  and  the  local 
contribution  to  DA  for  each  value  of  /3  is  calculated  using  the  expression 


XI  =  —  A 


77 


(18) 


where  the  factor  0.5  =  E3  (0).  The  factor  DA  is  then  given  by 

180° 


DVs=  SD  =  f  -  [°-5  -  ^Z)] 
0=0 


(19) 


11 


AEDC-TR-67-26 


These  DA's  are  next  summed,  as  previously  done  for  the  disk, 
to  give  the  man  quantities,  viz, 

a  =  £,  DA 

m  n 

with 


S  =  constant  -  O.ln 
St  =  constant  =  0.1m 

Once  these  are  found,  the  solution  for  y(Si)  is  given  by  the  matrix  inver¬ 
sion  technique  used  for  the  disk. 


SECTION  lit 

DISCUSSION  AND  CONCLUSIONS 


The  results  of  the  disk  calculation  are  shown  in  Fig.  1.  Also  shown 
for  comparison  are  the  flux  densities  for  free-molecule  flow  (R/A — «-0) 
and  continuum  flow  (R/A — ►  «>).  As  would  be  expected,  the  greater  shield¬ 
ing  occurs  at  the  axis  of  the  disk,  and  decreases  with  distance  from  the 
axis.  One  interesting  result  is  that  the  shielding  at  a  given  distance 
beyond  the  edge  of  the  disk  is  a  maximum  for  a  particular  value  of  R./A 
which  depends  upon  the  distance  from  the  disk. 

A  typical  result  of  the  cone  calculation  is  shown  in  Fig.  2  for 
R/A  =5.0,  The  validity  of  the  flow  model  used  for  the  calculation  is 
probably  questionable  for  this  relatively  high  value  of  R/A,  but  this 
particular  result  is  used  to  illustrate  the  main  features  of  the  calculation. 
The  most  important  result  is  that  the  point  of  the  cone  does  not  experience 
a  free-molecule  flux  density  except  in  the  case  of  a  vanishingly  small  cone 
angle.  Also,  except  for  the  region  adjacent  to  the  base  of  the  cone,  the 
incident  flux  level  is  identical  for  any  cone  length  in  terms  of  x/A.  This 
indicates  that  most  of  the  shielding  of  a  point  on  the  cone  is  attributable  to 
molecules  re-emitted  from  the  body  upstream  of  the  point.  Finally,  the 
lesser  shielding  of  a  point  by  all  locations  aft  of  that  point  is  related  to 
the  sudden  increase  in  flux  density  near  the  aft  end  of  the  cone.  Also 
shown  for  comparison  is  the  corresponding  result  for  the  normal  disk 
( 9C  -  90  deg). 

Figure  3  shows  the  effect  of  varying  R/A  for  a  particular  cone  angle. 
As  would  be  expected,  the  incident  flux  density  decreases  uniformly  as 
R/A  increases. 
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The  results  obtained  for  all  the  calculations  are  tabulated  in  Table  I. 
The  disk  is  taken  as  a  cone  of  0C  =  90  deg.  Table  II  gives  the  results 
for  the  influence  of  the  disk  in  the  plane  of  but  beyond  the  edge  of  the 
disk. 


It  is  of  interest  to  interpret  these  results  in  terms  of  quantities  use¬ 
ful  in  aerodynamics.  Since  in  free-molecule  flow,  pressure,  shear 
stress,  and  heat  transfer  are  proportional  to  free-stream  density,  the 
tabulated  quantity,  local  may  be  interpreted  in  those  terms,  viz, 


V  =  p.  p 


FM 


=  q  q 


FM 


KM 


The  drag  coefficient  ratio  is  given  by 


Cd  cDrM 


/ 


Normal  Area 


L  A"  i  dAn 


(20) 


which  is  also  equal  to  Q/QpM>  the  total  heating  rate  ratioed  to  the  free- 
molecular  value.  These  results  are  shown  in  Fig.  4  as  a  function  of  R/A. 
It  is  noted  that  as  R/A — -®,  C^/Cg^^ — >-0.  5,  which  is  approximately 
correct  for  a  blunt  body. 


Comparison  of  these  results  with  experimental  data  requires  that 
the  mean  free  path  of  the  re -emitted  molecules  be  defined.  In  Ref.  3  it 
is  shown  that 


A  =  A 


G 


1  +  \1 


where  C  =  2  (or  M„  k  6.0 


ay 

\f  9T7 


A. 

T, 


(21) 


Also,  Knudsen  number,  D  =  C y  2  ■  ReM 

where  C  is  a  constant  and  D  =  2R.  This  gives 


It  A 


lie 


2  Cy  M 
1  <> 


V  2 


(22) 


An  experimental  value  of  C  -  2.  446  is  given  in  Ref.  3  for  spheres  in  a 
hypersonic,  high-enthalpy  nitrogen  stream. 


It  is  of  interest  to  compare  the  present  results  with  similar  calcula¬ 
tions  by  Kogan  and  Degtyarev  published  in  Ref.  5.  Their  results  are 
based  on  an  approximate  solution  of  the  Boltzmann  equation  by  the  Monte 
Carlo  method,  and  allow  the  existence  of  a  molecular  boundary  layer  as 
Kogan  previously  postulated  in  Ref.  6.  A  comparison  between  their  cal¬ 
culations  and  the  present  results  is  shown  in  Fig.  5, 
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Fig.  1 


Fig.  5  Incident  Flux  Distributions  by  Two  Methods 
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TABLE  I 

TABULATED  VALUES  OF  INCIDENT  FLUX  DENSITY 


0 

0.1 

0.2 

0.3 

0.4  0.5 

R/X  -  0.05 

0.6 

0.7 

0.8 

0.9 

1.0 

5° 

0.9996 

0.9947 

0,9910 

0.9876 

0.9837 

0.9798 

0.9760 

0.9721 

0.9820 

0.9643 

0.9649 

10° 

0 . 9988 

0 . 9938 

0 . 9899 

0.9862 

0.9820 

0.9778 

0.9737 

0.9695 

0.9654 

0.9613 

0.9618 

20° 

0.9953 

0 . 9902 

0.9862 

0.9822 

0.9780 

0.97  38 

0.9698 

0.9656 

0.9619 

0.9583 

0.9594 

30° 

0.9891 

0.9840 

0.9801 

0.97  6  0 

0.9720 

0.9680 

0.9644 

0.9608 

0,9575 

0.9547 

0.9567 

90° 

0.9166 

0.9168 

0.9175 

0.9186 

0.9201 

R/X 

0.9222 

■=■  0.1 

0.9249 

0.92B4 

0.9328 

0.9387 

0.9476 

5° 

0.9993 

0.9897 

0.9824 

0.97  58 

0.9683 

0.9610 

0.9538 

0.9465 

0.9393 

0.9323 

0.9338 

10° 

0.9982 

0.9882 

0,9806 

0.9734 

0.9655 

0.9576 

0.9499 

0.9421 

0.9345 

0.9271 

0.9282 

20° 

0.9924 

0.9822 

0.9746 

0.9669 

0.9590 

0.9511 

0.9437 

0.9363 

0.9291 

0.9225 

0.9248 

30° 

0.9818 

0.9718 

0.9643 

0.9566 

0.9490 

0.9416 

0.9349 

0.9282 

0.9222 

0.9171 

0.9210 

90° 

0.8546 

0.8550 

0.8561 

0.8579 

0.9607 

R/X 

9.8643 

“  0.2 

0.8691 

0.8751 

0.8830 

0.8937 

0.9102 

5° 

0.9988 

0.9799 

0 . 9660 

0.9537 

0.9399 

0.9267 

0.9138 

0.9011 

0.8887 

0.8768 

0.8806 

10° 

0.9973 

0.9778 

0.9633 

0.9500 

0.9354 

0.9212 

0.9075 

0.8941 

0.8809 

0.8684 

0.8716 

20° 

0 . 9888 

0.9690 

0.9546 

0.9407 

0.9262 

0.9122 

0.8991 

0.8862 

0.8740 

0.B629 

0.8679 

30° 

0.9721 

0.9529 

0 . 9390 

0.9251 

0.9115 

0.8983 

0.8867 

0.87  52 

0.8648 

0.8564 

0.8641 

90° 

0.7608 

0.7673 

0.7690 

0.7719 

0.7761 

R/X 

0 '.  7  8 1 9 
-  0.5 

0.7894 

0.7992 

0.8121 

0.8301 

0.8594 

5° 

0.9974 

0.9525 

0.9218 

0.8962 

0.9686 

0.8436 

0.8202 

0.7983 

0.7778 

0.7590 

0.7718 

10° 

0.9956 

0.9494 

0.9177 

0.8902 

0,8611 

0.8345 

0 . 8098 

0.7866 

0.7649 

0.7450 

0.756  5 

20° 

0 . 9839 

0.9374 

0.9067 

0.8785 

0.8503 

0.8244 

0.8011 

0.7791 

0.7588 

0.7412 

0.7555 

30° 

0.9586 

0.9144 

0.8855 

0.8582 

0.8321 

0.8080 

0.7B74 

0.7678 

0.7506 

0.7372 

0.7559 

90° 

0.6330 

0 . 6337 

0.6360 

0.6398 

0.6455 

R/X 

0.6535 

-  1.0 

0,6642 

0.6787 

0.6990 

0.7292 

0.7866 

5° 

0.9952 

0.9119 

0.8613 

0.8227 

0.7830 

0.7496 

0.7203 

0.6943 

0.6710 

0.6507 

0.6789 

10° 

0.9933 

0.9077 

0.8559 

0.8144 

0.77  28 

0.7373 

0.7065 

0.6789 

0.6543 

0,6328 

0.6587 

20° 

0.9800 

0.8946 

0.8453 

0.8033 

0.7638 

0.7298 

0.7009 

0.6749 

0.6518 

0.6325 

0.6615 

30° 

0.949B 

0.8699 

0.8244 

0.7845 

0.7484 

0.7171 

0.6915 

0.6683 

0.6484 

0.6335 

0.6675 

90° 

0 . 5526 

0.5532 

0.5551 

0.5585 

0.56  36 

R/X 

0 . 57 10 

-  2.0 

0.5814 

0 . 5964 

0.6191 

0.6561 

0.7451 

5° 

0.9914 

0.8456 

0.7745 

0.7273 

0.6816 

0.6473 

0.6194 

0.5957 

0.57  54 

0.5583 

0.6099 

10° 

0.9893 

0.8398 

0.7675 

0.7166 

0.6690 

0.6323 

0.6027 

0.5774 

0.5556 

0.5371 

0.5850 

20° 

0.9746 

0.8268 

0.7595 

0.7086 

0.6641 

0.0295 

0.6017 

0.5777 

0 . 5570 

0.5400 

0 . 5906 

30° 

0.9409 

0.8045 

0.7436 

0.6955 

0.6556 

0.6239 

0. 5992 

0.5775 

0.5591 

0.5456 

0.6011 

90° 

0.5114 

0.5117 

0.5126 

0.5144 

0.5172 

R/X 

0 . 5217 

=  5.0 

0.5285 

0 , 5393 

0.5580 

0 . 5931 

0.7277 

5° 

0.9820 

0.7189 

C  .  6464 

0.6090 

0 . 57  29 

0.5495 

0.5310 

0.5147 

0.5001 

0.4871 

0 . 57  54 

10° 

0.9809 

0.7108 

0,6379 

0.5966 

0.5584 

0.5325 

0.5121 

0.4940 

0,4776 

0.4631 

0,5447 

20° 

0.9651 

0.7008 

0.6353 

0 . 5933 

0.5584 

0.5341 

0.5147 

0.4974 

0.4816 

0.4677 

0.5512 

30° 

0.9269 

0.6868 

0 . 6287 

0.58  9  0 

0.5584 

0.5364 

0.5188 

0.5027 

0 . 4880 

0.4758 

0.5634 

90° 

0.5002 

0.5003 

0 . 5003 

0.5005 

0 . 5009 

0. 5015 

0.5029 

0.5057 

0 . 5130 

0 . 5246 

0.7125 
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TABLE  II 

INCIDENT  FLUX  DENSITY  IN  PLANE  OF  DISK 


1.1 

1.2 

1.3 

1.4  1,5  1.6 

R/X  =0.05 

1.7 

1.8 

1.9 

2.0 

0.9544 

0.9582 

0.9611 

0.9636 

0.9656  0.9674 

R/X  =0.1 

0.9689 

0.9702 

0.9715 

0.9725 

0.9235 

0.9307 

0.9364 

0,9410 

0.9449  0.9483 

R/X  =0.2 

0.9512 

0.9538 

0.9561 

0.9581 

0.8845 

0.8977 

0.9080 

0.9163 

0.9233  0.9292 

R/X  =0.5 

0.9343 

0.9388 

0.9438 

0,9463 

0.8412 

0.8674 

0.8869 

0.9021 

0.9144  0.9346 

R/X  =1.0 

0.9331 

0.9403 

0.9465 

0.9518 

0.8370 

0 . 87  56 

0.9022 

0.9216 

0.9362  0.9476 

R/X  =2.0 

0.9565 

0,9637 

0.9694 

0.9742 

0.86  53 

0.9133 

0.9415 

0.9594 

0.9712  0,9792 

R/X  =  5,0 

0.9848 

0 . 9888 

0.9917 

0.9938 

0.9294 

0,9712 

0.9871 

0.9939 

0.9970  0,9985 

0,9992 

0.9996 

0.9998 

0.9999 
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13  ABSTRACT 

The  calculations  presented  herein  are  based  on  a  first-collision 
model  which  allows  collisions  between  the  free-stream  molecules  and 
the  molecules  which  are  re-emitted  from  the  surface  of  the  body.  The 
net  effect  of  these  collisions  is  to  partially  shield  the  body  from 
the  free  stream,  reducing  both  the  drag  and  heat-transf er  coefficients 
from  the  corresponding  values  experienced  in  free-molecule  flow. 
Free-stream  Mach  number  is  taken  to  be  essentially  infinite  and  the 
molecules  are  assumed  to  be  re-emitted  from  the  body  surface  at  the 
most  probable  velocity,  instead  of  possessing  a  velocity  distribution. 
These  assumptions  enable  the  distribution  of  flux  to  a  given  body  to 
be  expressed  as  a  function  of  the  degree  of  rarefaction,  as  repre¬ 
sented  by  the  appropriate  Knudsen  number.  This  Knudsen  number  is 
composed  of  a  characteristic  body  dimension  and  the  mean  free  path  of 
the  re-emitted  molecules  relative  to  free-stream  molecules.  The 
integral  equation  expressing  the  incident  flux  distribution  on  a 
general  body  is  developed,  and  solutions  are  presented  for  the  disk 
normal  to  the  free  stream  and  for  sharp  cones  of  various  apex  angles, 
at  zero  angle  of  attack.  These  flux  distributions  are  then  integrated 
to  give  drag  coefficients  ratioed  to  the  corresponding  free-molecule 
values. 
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